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We report on a program to compute the electromagnetic form factors of mesons. We discuss the 
techniques used to compute the pion form factor and present results computed with domain wall 
valence fermions on MILC asqtad lattices, as well as with Wilson fermions on quenched lattices. 
The methods can easily be extended to p — > 77r transition form factors. 
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Q , I. INTRODUCTION 

The pion electromagnetic form factor is often considered a good observable for studying the onset, with increasing 
energy, of the perturbative QCD regime for exclusive processes. It is believed that, because the pion is the lightest 
and simplest hadron, a perturbative description will be valid at lower energy scales than predictions for heavier and 
more complicated hadrons such as the nucleon 

A pseudoscalar particle has only a single electromagnetic form factor, F(Q 2 ), where Q 2 is the four-momentum 
OO . transfer, and furthermore at Q 2 = 0, this form factor is normalized to the electric charge of the particle, F(Q 2 = 0) = 1; 
the magnetic form factor vanishes. Thus in this paper, we will be measuring the form factor of the positively charged 
7T + . The experimentally observed behavior of the form factor at small momentum transfer is well described by the 
vector meson dominance (VMD) hypothesis S LI 

O F^iQ 2 ) ~ 2 for Q 2 <^m 2 VMD . (1) 
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The current experimental situation is presented in Fig. ^ 

For Q 2 < 0.28 GeV 2 , Amendolia et al. [BJ have determined the form factor to high precision from scattering of very 
| high energy pions from atomic electrons in a fixed target. For higher Q 2 , the form factor has been determined from 
quasi-clastic scattering from a virtual pion in the proton In this case, the extracted values for the form factor 

must depend upon some theoretical model 0, UnL for extrapolating the observed scattering from virtual pions to 
. ,— i , the expected scattering from on-shell pions. As the models have become more sophisticated, the earlier data |(| Q 
■ have been reanalyzed [l^ for consistency. Shown in Fig. are the results of some model calculations which seem to 
H \ cover the range of existing predictions [l^, HJ, Hal ■ 

Given the dominance of the rho meson resonance in the time-like region (Q 2 < 0) of the pion form factor, perhaps 
it is not too surprising that the space-like form factor is well described at low Q 2 > by a VMD-inspired monopole 
form with only a contribution from the lightest vector resonance (mvMD ~ Tip ~ 0.77 GeV). What is striking is that 
it accurately describes all experimental data even up to scales of Q 2 > 1 GeV 2 . Furthermore, VMD predicts that the 
form factor should scale as F 7T (Q 2 ) ~ 1/Q 2 for Q 2 3> m 2 /MD , the same scaling predicted at asymptotically high Q 2 in 
perturbative QCD [l7L lls[ . One crucial fact which makes the pion form factor an ideal observable for studying the 
interplay between perturbative and n on-p ert urbative QCD is that its asymptotic normalization can be determined 



from pion decay p.a ]M UA 121 XM IM U. 
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FIG. 1: Summary of experimental data for the pion electromagnetic form factor; shaded regions indicate expected sensitivities 
and coverage of future results. The lines indicate theoretical calculations of the form factor, as described in the text. 

Higher order perturbative calculations of the hard contribution to the form factor [26l I2H [28l |29| do not vary 
significantly from this value, and are shown in Fig. ^ At the largest energy scale where reliable experimental 
measurements have so far been obtained, around Q 2 ~ 2 GeV 2 , the data are 100% larger than this pQCD asymptotic 
prediction. 

This situation raises many questions. At what scale does the form factor vary with Q 2 as predicted by pQCD? 
However, merely observing the proper Q 2 dependence is insufficient as Eq. (JJJ has the same asymptotic Q 2 dependence 
as pQCD but is numerically about twice as large. How rapidly will the data approach the pQCD prediction and at 
what scale will pQCD finally agree with the data? These are questions which Lattice QCD calculations are ideally 
suited to address, provided we can get reliable results for momentum transfer on the order of a few to several GeV 2 . 

Early lattice calculations validated the vector meson dominance hypothesis at low Q 2 |30ll3l| . Recent lattice results 
[32T l33l l3~iL I35L l3?l I37J , including some of our own preliminary results [3^, |3?| , have somewhat extended the range of 
momentum transfer, up to 2 GeV 2 , and the results remain consistent with VMD and the experimental data. 

II. LATTICE COMPUTATION OF F„{Q 2 ) 

The electromagnetic form factor is obtained in lattice QCD simulations by placing a charged pion creation operator 
at Euclidean time ij, a charged pion annihilation operator at tf and a vector current insertion at t as shown in 
Fig. [3 A standard quark propagator calculation provides the two propagator lines that originate from U. The 
remaining quark propagator, originating from tf is obtained via the sequential source method: (1) completely specify 
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the quantum numbers, including momentum p^. of the annihilation operator to be placed at tf and (2) contract the 
propagator from t, to tf to the annihilation operator and use that product as the source vector of a second, sequential 
propagator inversion. The resulting sequential propagator appears as the thick line in Fig. [3 extending from U to t 
via tf. Given these two propagators, the diagram can be computed for all possible values of insertion position t and 
insertion momenta q; the initial momentum p ; is determined by momentum conservation = — q. 

Furthermore, with the same set of propagators, any current can be inserted at t and any meson creation operator 
can be contracted at i,. So, the diagram relevant to determining the form factor for the transition p + — » 7*7r + can 
be computed without further quark propagator calculations. By applying the sequential source method at the sink, 
the trade-off is that the entire set of sequential propagators must be recomputed each time new quantum numbers 
are needed at the sink, particularly p^. 




FIG. 2: The quark propagators used to compute the pion form factor. 



We can extract the pion energies E„ (p) using standard lattice techniques of fitting pion correlation functions from 
which we can compute the momentum transfer 



-Q 2 =[E v { V f)-E«{ Vi )] 2 -{Vf-vtf 



(3) 



which should be non-positive if the pion spectral function is well-behaved. Since the largest Q 2 for a given |p| 2 occurs 
in the Breit frame, P/ = — P,, it is important to choose a non-zero P/ to achieve large momentum transfer; indeed 
Pf = yields a vanishing Q 2 for all q in the chiral limit. 
The form factor, F(Q 2 ), is defined by 



<7T( P/ ) |^(0)| 7T( Pi )) 



continuum 



Zy <7T( P/ ) |^ at (0)| 7T( Pi )> = F{Q 2 ){ Pl +p f ) fi 



(4) 



where V^ at (a;) is a particular lattice discretization of the continuum vector current. In this work, local and point-split 



vector currents were used: 



(5) 
(6) 



Zy is the corresponding matching factor between the lattice current and continuum; for Wilson fermions, the point- 
split current is conserved on the lattice: Z v ' s ' = 1 (Wilson). The three-point correlation function depicted in Fig.0 
is given by 



T^(U,t,tf,p t ,p f )= £ e -( x /- x )-P/ e -*(x-xO- P< Mxf)V ^ A{Xi) 



(7) 



where 4>\{x) and 4>b{%) are creation and annihilation operators with pion quantum numbers. The A and B indicate 
that different operators may be used at the source and sink, i. e. smeared source and point sink or pscudoscalar source 
and axial vector sink. 

Inserting complete sets of hadron states, and requiring U t <C tf, gives 



-(t f -t)E„(p 



r^(< i ,t,t / ,p i , P/ ) -» (o|<mo)Wp/)) 

Similarly for the two-point correlator, with fcj < t/, 



<7r( P/ ) |^(0)| 7r(p 4 )> 



-(t-t i )E 7 ,(p i 



TT(Pi) 4(0) 



r£*(u,t f , P ) -> (o |^(o)| tt(p)) e " ( lC'^ (P) ( 



(8) 



(9) 
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We use both the pseudoscalar density, (jy^{x) = tp(x)j5%jj(x), and the temporal component of the axial vector 
current, <fi( 2 \x) = ^(x^^j^i/j^x), as pion interpolating operators, constructed from both local and smeared quark 
fields, denoted by L and S respectively. We adopt gauge-invariant Gaussian smearing 



b(x) 



1 + -VU) b(x) 



N 



(10) 



where uj and N are the tunable parameters used to specify the smearing radius; the flavor structure is suppressed. The 
respective merits of these interpolating operators will be discussed later. For these operators, wc define the following 
amplitudes 



0^(0) ^p: 
4 2) (o)|^(p: 

4 2) (0) Ap] 



= Z 



Z ( s'(\p\) 



= z 



(2), 



(2) 

The overlap of the operator <p L has trivial p dependence arising from the Lorentz structure of the operator 



(11) 
(12) 

(13) 
(14) 

(15) 



However, the introduction of an additional three-dimensional scale introduces non-trivial p dependence for the smeared 
overlaps Zg and Zg . 

We employ two methods to determine the form factor F 7T {Q 2 ). The first method, which wc call the fitting method, 
involves a fit of the relevant two- and three-point functions to simultaneously extract the form factor, the energies 
E^(p) and the amplitudes Z(p) in a single covariant, jackknifcd fit. 

The second method, which we call the ratio method, starts by determining the energies £V(p), either by fits to the 
correlators at non-zero momentum or from a dispersion relation, and then constructing the following ratio which is 
independent of Zr ', Zs(|p|) and all Euclidean time exponentials for sufficiently large temporal separations: 



F{Q\t) 



r^(t ( ,t ) t /j p < ,p / )r^(t i) t,p / ) / 2Z v E n (p f ) 



T% ( ti ,t, Pi )TCB (U,t f ,p f ) \E„ ( Pi ) + E w (p f ) 



(16) 



where the indices A, B and C can be either L (local) or S (smeared). As part of our program, we expect to determine 
the relative merits of each extraction method. 



III. SIMULATION DETAILS 



Our first calculations were performed on quenched configurations generated with the Wilson gauge action at (3 = 6.0, 
corresponding to a -1 « 2 GeV. The propagators were computed using the unimproved Wilson fermion action with 
Dirichlet boundary conditions in the temporal direction, and periodic boundary conditions in the spatial directions. 
For Wilson fermions at these lattice spacings, the exceptional configuration problem is rather mild particularly when 
compared to that for the non-perturbatively improved clover action. This enabled us to reach pion masses of 300 McV 
without observing any exceptional configurations, whereas computations using the clover action are limited to pion 
masses in excess of roughly 400 MeV [^3 due to the systematic errors caused by the frequent occurrence of exceptional 
configurations. A detailed listing of the simulation parameters is provided in Tab. [J 

The pion masses attainable in quenched domain-wall fermion (DWF) calculations are limited only by finite volume 
effects and available computing power. So far, however, quenched DWF computations of the pion form factor have 
only explored pion masses down to 390 MeV, at a" 1 ~ 1.3 GeV [3j|. An important advantage of DWF fermions is 
that they are automatically 0(a)-improved. An alternative approach using twisted-mass QCD (tmQCD), in which 
Wilson fermion results are also 0{a) improved with just double the effort^^i h as been explored separately [34j. 

The Lattice Hadron Physics Collaboration (LHPC) has been performing unquenched hadron structure calculations 
using MILC Nf = 2 + 1 and Nf = 3 configurations generated with staggered asqtad sea quarks |4ll l42l \I$L E3. HEj. 
Valence propagators were computed using domain wall quarks with domain wall height m = 1.7 and extent L s = 16 
of the extra dimension and Dirichlet boundary conditions imposed 32 timeslices apart. The MILC configurations were 
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TABLE I: Simulation details for quenched Wilson fermion calculations at a 1 ~ 2 GeV 



k volume 


am p am-r: m^/m p m w L 




0.1480 16 3 x 32 
U.loUU 10 X oZ 

0.1520 16 3 x 32 
0.1540 16 3 x 32 


0.7187(39) 0.6752(45) 0.943(10) 10.80(7) 

U.DOO<(0/J U.0oo4^4oj U.yULUloj i).OO(l) 

0.5540(83) 0.4851(67) 0.876(14) 7.76(11) 
0.4682(124) 0.3752(73) 0.801(23) 6.00(11) 


0.82725(51) 

U. i o4eo^40 ) 

0.74197(54) 
0.70459(39) 


0.1555 24 3 x 32 
0.1563 24 3 x 32 


0.4209(88) 0.2613(29) 0.621(14) 6.27(7) 
0.4014(68) 0.1921(29) 0.479(10) 4.61(7) 


0.65676(43) 


0.1566 32 3 x 48 


0.3724(145) 0.1629(36) 0.437(19) 5.21(12) 


0.65553(14) 



TABLE II: Simulation details for domain wall fermion calculations on 20 3 x64 dynamical MILC asqtad lattices at a 1 
1.6 GeV. 



am ud 


am s 


arrival 


m p (MeV) 


TUn (MeV) 


m w /m p m^L 


rjXoc ^p.S. 


0.01 
0.05 
0.05 


0.05 
0.05 
0.05 


0.01 
0.05 
0.081 


956(22) 
955(19) 
1060(14) 


318(3) 
602(5) 
758(5) 


0.333(8) 3.97(4) 
0.630(12) 7.68(6) 
0.715(10) 9.66(6) 


1.0714(55) 1.1098(23) 
1.0890(55) 1.0835(37) 
1.1199(14) 1.0833(13) 



HYP blocked |4g| before the valence propagators were computed to avoid unacceptably large residual chiral symmetry 
breaking. For our pion form factor calculation we chose configurations separated by at least 12 HMD trajectories. 
In some of the other LHPC calculations, if measurements were computed on configurations separated by 6 HMD 
trajectories, opposite halves of the lattice were used to minimize autocorrelations. The simulation parameters for the 
DWF computation are provided in Tab. |H| A detailed study of the physi cal p roperties of light hadrons composed of 
staggered quarks computed on these lattices has recently been completed [47j . 



IV. RESULTS 
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FIG. 3: Results for the pion form factor as a function of Q 2 for each of the available k values, together with experimental 
determinations 0. The curve shows the vector-meson-dominance prediction, using the experimental determination of the p 
mass. 

Results for the quenched Wilson form factor computed by the fitting method at each value of the quark mass are 
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FIG. 4: Pion (left) and rho meson (right) dispersion relation vs. continuum (upper) and lattice (lower) expectations curves. 

shown in Fig. |2| Also indicated are experimental data points and the vector meson dominance (VMD) prediction 
using the experimental value for the p meson mass. While the data tend in the correct direction with decreasing 
pion mass, the reader may notice that the form factor for 300 MeV pions already lies below the physical curve. This 
suggests that the pole mass that would be obtained from a VMD fit to the lattice data would be considerably lower 
than the physical p meson mass. This is perhaps expected since it is known that 0(a) scaling violations in the vector 
meson mass computed with Wilson fermions yield an underestimate of the mass of roughly 20% 01 j the same amount 
needed to move the form factor points so as to lie above the continuum curve. 

Because we would like to compute the pion form factor at large momentum transfer, we have spent a substantial 
amount of effort on our domain wall data set in extracting the pion energies and amplitudes at relatively large 
momenta; the largest attainable momenta is constrained by the fineness of the lattice spacing. In the continuum 
limit, the pion dispersion relation should approach the continuum one 



K(0) 5 



(17) 



At a non-zero lattice spacing, a study of free lattice bosons suggests the lattice dispersion relation 



E(pf 



(18) 



where E and p are the "lattice" energy and momentum respectively: 



E = 2 sinh 



E 



2 sin 



(19) 



Eqs. I|17J) and l|18(l agree in the small- momentum limit. 

In Fig. 01 we show the measured "lattice" energies against the lattice momenta at each of our quark masses, together 
with curves representing the continuum and lattice dispersion relations. We see that both dispersion relations provide 
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FIG. 5: Results for the local (L) and smeared (S) pseudoscalar (A, A) and axial vector (■, □) pion source amplitudes (p) 
used in this study. 



a reasonable representation of the data, although there may be a slight flattening of the data against the continuum 
curve at higher momenta. These results suggest that directly fitting all the data to either dispersion relation, thereby 
reducing the number of fit parameters needed to extract the form factor, may improve the relative signal to noise of 
the remaining parameters. This may help dramatically in the ratio method, where the only fit parameters are the 
form factor and the energies. 

In the fitting method, one must reliably extract not only the energies, but the amplitudes, at high momenta. In 
Fig. [S] we present the four amplitudes of Eqs. Ijllfl through (|14H that we estimate from the four two-point correlators 
we measure: smeared-smeared and smeared-local for both pseudoscalar-pseudoscalar and axial-axial operators. In the 
fitting procedure, all four correlators are constrained to have the same energy. From the figure, we can see that our 
expectations of Zjp oc const and zf \p) oc -E^p) are consistent with the data. We can also see from Z^\p = 0) 
that the smeared pseudoscalar operator has a strong overlap with the zero momentum pion relative to the axial- vector 
operator. 

We plot in Fig.[S]just the relative error estimates for the amplitudes. This clearly demonstrates that the statistical 
noise inherent in a particular source or sink operator need not be correlated with the magnitude of the amplitude. 
In particular, we see that beyond the few smallest momenta, the pseudoscalar and axial-vector amplitudes are of the 
same order, but the inherent noise of the pseudoscalar operator is unacceptably large for higher momenta. 

Using the lattice dispersion relation and the ratio method described above and in our previous work, we have 
computed in Fig.[fJlthe pion form factor for two of the dynamical pion masses. Fits of the data to the monopole form 
of Eq. arc also shown, where shaded regions correspond to jackknife error bands and the central values for the 
pole masses are given in the legend. The data in Fig. [B] were computed with pseudoscalar pion sink operator fixed at 
momentum = (1,0,0). Tab. Illll summarizes the fits for all available combinations of m V ai and p^-. Unconstrained 
linear extrapolation of the two most reliable data points, m va i = 0.01 and 0.081 with p^ = (1, 1,0), gives an estimate 
of mvMD = 0.868(65) GeV in the zero quark mass limit. From this we estimate the mean square charge radius of the 
pion to be ( r2 ) VMD = 0.310(46) fm , which is significantly below the experimental value of {r 2 ) cx t = 0.439(8) fm 2 . 




12 3 



d (GeV) 2 

FIG. 6: Pion electromagnetic form factor for fixed sink momentum p^- = (1,0,0) computed by the ratio method l38l l3Sl| and 
imposing the lattice dispersion relation, Eq. (1181 . Shaded regions are jackknife error bands for VMD fit. 



TABLE III: Results of fit of form factor data to the VMD monopole ansatz of Eq. Q. 



am va i 


P/ 


ijjvmd (GeV) 


X 2 /d.o.f. 


d.o.f. 


■^Vi confs 


0.01 


(0,0,0) 








51 




(1,0,0) 


0.888(56) 


1.37(72) 


13 


251 




(1,1,0) 


0.30(20) 


1.7(1.2) 


11 


106 


0.05 


(0,0,0) 


1.278(87) 


2.8(1.3) 


20 


104 


0.081 


(0,0,0) 


1.192(93) 


3(17) 


15 


49 




(1,0,0) 


1.030(73) 


2.1(1.7) 


22 


70 




(1,1,0) 


1.022(87) 


3.6(2.2) 


23 


73 



More dynamical quark masses, particularly lighter ones, will be needed to perform a proper chiral extrapolation. 

We were unable to compute statistically significant form factors with the same pseudoscalar sink operator for 
"ivai = 0.01 and = (1, 1, 0) by the ratio method, apparently due to the poor signal-to-noise inherent in the overlap 
of the chosen sink operator with the higher momentum state. For m va i = 0.01 and p* = (0, 0, 0) the signal was rather 
better but 51 sequential propagators were too few to allow for a stable fit to the VMD ansatz. We expect that an 
axial vector sink operator would be a better choice for future calculations at non-zero sink momentum. Some of us are 
also extending our group-theoretical construction of extended baryon operators 0, IBTl l52l l5^ | to include mesonic 
quantum numbers. 
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V. CONCLUSIONS 

From our quenched Wilson form factor results, we find that both the ratio method and the fitting method are useful 
tools for computing the pion form factor. Each method has different systematic errors, so the extent to which both 
agree should give confidence that the systematic errors are small and well understood. A comparison of the Wilson 
form factor results with the experimental pion form factor data suggests that the former have large, and expected, 
discretization errors. 

Domain- wall fcrmions are free of 0(a) discretization uncertainties to the extent that any residual mass is eliminated. 
Thus our dynamical fermion results are obtained using domain-wall fermions for the valence quarks computed on 
asqtad lattices. From a careful analysis of the pion spectrum, we recognize the importance of using a large basis of 
pion operators so that we can identify at least one operator at each pion momentum whose overlap with the pion has 
a reasonable signal-to-noise ratio. The local axial vector operator appears to be a better choice than the pseudoscalar 
operator for fitting higher momentum states since its overlap with a given state increases in proportion with the pion 
energy and without any degradation of the signal-to-noisc ratio. Our analysis enables us to obtain the pion form 
factor for a pion mass approaching 300 MeV at a range of Q 2 commensurate with the current experimental program. 

Given our existing analysis framework and the propagators and sequential propagators that we have already com- 
puted, we can also compute the p + — » 771-+ transition form factor. We hope to complete this analysis in the near 
future. 
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